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QQ modes in the Quark-Gluon Plasma 
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Abstract. We study the evolution of heavy quarkonium states with temperature in a Quark-Gluon Plasma 
by evaluating an in-medium QQ T— matrix within a reduced Bethe-Salpeter equation in S— and P— wave 
channels. The interaction kernel is extracted from finite-temperature QCD lattice calculations of the singlet 
free energy of a QQ pair. Quarkonium bound states are found to gradually move across the QQ threshold 
after which they rapidly dissolve in the hot system. We calculate Euclidean-time correlation functions and 
compare to results from lattice QCD. We also study finite- width effects in the heavy-quark propagators. 

PACS. 25.75.Dw - 12.38.Gc - 24.85. -hp - 25.75.Nq 



1 Introduction 

Bound states of heavy (charm and bottom) quarks {Q= 
b, c) are valuable spectroscopic objects in Quantum Chro- 
modynamics (QCD) When embedded into hot and/or 
dense matter a large class of medium modifications can be 
studied, including (Debye-) color-screening of the QQ in- 
teraction, dissociation reactions induced by constituents 
of the medium, and the change in thresholds caused by 
mass (or width) modifications of open heavy-fiavor states. 
Lattice QCD (IQCD) calculations have made substantial 
progress in characterizing in-medium quarkonium proper- 
ties from first principles. In particular, it has been found 
that ground state charmonia |2||3M] and bottomonia |5j do 
not dissolve until significantly above the critical tempera- 
ture, Tc, which has been supported in model calculations 
based on potentials extracted from IQCD, using either a 
Schrodinger equation for the bound-state problem ,fe7n8 , 
, or a T- matrix approach which additionally accounts for 
scattering states . More reliable comparisons to IQCD 
can be performed using (spacelike) Euclidean-time cor- 
relation functions |11II12| . which are readily evaluated in 
IQCD. The conversion of (timelike) model spectral func- 
tions requires a description not only of the bound state 
part of the spectrum but also its continuum and thresh- 
old properties. 

In the present work we evaluate Euclidean correlation 
functions for charmonium and bottomonium in a T— matrix 
approach. The basic input consists of in-medium QQ po- 
tentials extracted from IQCD, inserted in a scattering equa- 
tion to calculate the in-medium QQ T— matrix This 
incorporates bound and scattering states on an equal foot- 
ing (based on the same interaction), and additionally en- 
ables a straightforward implementation of in-medium single- 
particle (quark) properties via self-energy insertions in the 
two-particle Green's function. 



2 Scattering equation and bound states 

The T-matrix equation for QQ scattering in the center of 
mass frame and in partial wave basis reads jlO) . 

TiiE;q\q)^Viiq',q) 

- - / dkk^VM' ,k)GQQ{E-k)Ti{E-k,q) , 
Jo 

(1) 

which follows from a standard 3-dimensional reduction of 
the Bethe-Salpeter equation ^qq{^t^) denotes the 
intermediate two-particle propagator including quark self- 
energies (S) and Pauli blocking. The T— matrix equation 
(Q) is solved with the Haftel-Tabakin algorithm in 
which the integral equation is solved by discretizing 3- 
momentum, "^^^i ^{E)ikT{E)kj — Vij, and subsequent 
matrix inversion. The zeroes of det ^{E) for E < Eth de- 
termine heavy quark-antiquark bound states. 

The quark selfenergy, S, encodes the interactions with 
(light) quarks and gluons from the heat bath [THj. In the 
present work we consider a fixed heavy-quark mass mq 
(i.e.. Re — 0) together with a small imaginary part, 
ImS = —0.01 GeV, for numerical purposes. Effects of 
temperature dependent heavy-quark masses and widths 
are investigated in Ref. [T5l . 

The kernel of the scattering equation, V, can be es- 
timated from the IQCD heavy-quark free energies, Fi, 
even though there is still an ongoing discussion on how 
to properly do it. Here we identify the heavy-quark po- 
tential with the color-singlet internal energy, Ui = Fi — 
T , which reproduces ground-state charmonium disso- 
ciation temperatures as found in lattice analysis of spec- 
tral functions [7||Sl|^10j. In Fig. ^ we show the QQ poten- 
tial, V{r, T) = Ui{r, T) -Ui{r ^ oo, T), as obtained from 
a fit to the IQCD color-singlet free energy data ^HIj pre- 
viously employed in JOI- The potential evolves smoothly 



2 



D. Cabrera, R. Rapp: QQ modes in the Quark-Gluon Plasma 




r[fm 



Fig. 1. QQ potential for several temperatures above Tc based 
on the color-singlet internal energy. 
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Fig. 2. T— matrix for S— wave cc scattering based on the 
potential in Fig. Also shown is det^ (dashed line, arbi- 
trary units). From left to right (up-down) the temperatures 
are (1.1, 1.5, 2.0, 2.5, 3.0, 3.3) T^. 



with temperature, decreasing both in magnitude and range. 
Different parameterizations of the lattice data imply siz- 
able uncertainties in the potential through the thermal 
derivative of the free energy. These uncertainties are stud- 
ied in Ref. jl5| by alternatively obtaining the heavy-quark 
potential from a direct fit to the IQCD internal energy 
data of Ref. ^iW t1) follows from a Fourier trans- 

form and partial-wave expansion. 



3 Quarkonium T— matrices in the QGP 

We start the calculation of the in-medium QQ T— matrices 
by fixing the heavy-quark masses so that the correspond- 
ing quarkonium ground states approximately agree with 
their vacuum masses for the lowest considered tempera- 
ture (T l.lTc), i.e., mc = 1.7 GeV and = 5.1 GeV. 
Fig. 121 summarizes the on-shell S*— wave cc scattering am- 
plitude as a function of the CM energy, for several tem- 
peratures. We do not include the hyperfine (spin-spin) in- 
teractions and therefore rjc (Vb) and J/tp (T) states are 
degenerate. At the lowest temperature, we recover the 
charmonium ground state at E Ri 3.10 GeV, whereas the 
first excited state has just about melted. As tem- 
perature increases, the J/ip{lS) gradually moves toward 
threshold and the T— matrix is appreciably reduced. The 
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Fig. 3. Same as in Fig. |21 but for 5*— wave bh scattering. 
From left to right and up to down the temperatures are 
(1.1,1.5, 1.8, 2.1, 2.7, 3.5) 

Table 1. Masses and binding energies (in [GeV]) for P— wave 
quarkonia as obtained from detT{E) — 0. 
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bound state survives up to ~ 3Tc, where it crosses the cc 
threshold and rapidly melts in the hot system. 

The S*— wave bb T— matrix exhibits two bound states 
at the lowest temperature {E 9.45, 9.95 GeV for r{lS), 
rjii and T(2S), -q'^^, respectively) and the remnant of a third 
one, which is (almost) melted in the medium, cf. Fig. |31 
The T{2S) moves across the bb threshold at T w 2.1 T^, 
whereas the IS' state survives in the QGP until much 
higher temperatures, beyond T « 3.5 Tc. 

We only find one P— wave bound state for the charm 
system at the lowest temperature, at i? « 3.4 GeV, which 
we associate with the Xc- The P— wave bb system exhibits 
two bound states at 1.1 Tc, with energies E = 9.90, 10.15 GeV, 
in good agreement with the nominal values for Xbi^P) and 
Xb(2P) in the vacuum. The latter moves beyond threshold 
at T « 1.3 Tc, and the (IP) state at T w 2.3 Tc. Masses 
and binding energies, {Eb — Eth — M), of the P— wave 
states are summarized in Tab.^for several temperatures. 



4 Spectral functions and Euclidean correlators 

The Euclidean-time correlation function is defined as the 
thermal two-point mesonic correlation function in a mixed 
Euclidean-time- momentum representation (here p = 0) , 



G(t,T) 



duj (t(uj, T) 



cosh[u;(r - /3/2)] 
sinh(w/3/2) 



(2) 
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where a is the spectral function obtained from closing ex- 
ternal legs of the in-mcdiuni T— matrix, schematically as 

G{E)= J Gqq+J GqqTGqq, aiE) ^ImGiE) . (3) 

As expected, the S"— wave charmonium spectral func- 
tion (left panel of Fig. Q reflects the bound states of 
the T— matrix, but the (non-perturbative) cc rescatter- 
ing also generates appreciable strength above threshold, 
where the remnant of first excited state (-0(25*)) can be 
still inferred at T = l.lTc. The corresponding correlation 
function (right panel of Fig. 0)) is normalized to a "re- 
constructed" correlator represented by a zero-temperature 
spectral function consisting of a ^-function like bound 
state spectrum and perturbative continuum with onset 
at Econt = 4.5 GeV IH]. The temperature evolution of 
the correlator is a combined result of a decrease in bind- 
ing energy of the bound states and the contribution of the 
non-perturbative continuum. The sizable drop at large r is 
in qualitative agreement with IQCD . The latter exhibit 
somewhat less reduction, leaving room for the effects of a 
threshold reduction with temperature in our T— matrix. 

The P— wave charmonium spectral function (Fig. 
left) exhibits one bound state (xc) just below the thresh- 
old, which rapidly melts into the continuum accompanied 
by a sizable threshold enhancement. The normalized cor- 
relator steeply rises in the low-r regime, due to: (i) non- 
perturbative rescattering and (ii) a larger threshold in the 
"reconstructed" correlator. While this is qualitatively con- 
sistent with IQCD, the evolution with temperature is op- 
posite, which could improve by a shift of strength to lower 
energies due to a decreasing heavy-quark mass with tem- 
perature. The bottomonium correlation functions follow a 
similar pattern as for the charmonium system. 

It turns out that the r-dependence of the (normal- 
ized) Euclidean correlators is rather sensitive to the "re- 
constructed" correlator used for normalization. However, 
5*— wave charmonium correlators from IQCD show rather 
little temperature dependence below ^l.STc. Thus, to re- 
duce ambiguities induced by the reconstructed correlator, 
we have normalized our correlators from the T— matrix to 
the result at T = 1.1 Tc calculated in the same approach, 
cf. Fig. IHI The r-dependence is now substantially weaker 
than in the right panel of Fig. 0] (limited to less than 10%), 
and the agreement with IQCD 3 is much improved. 

We have also studied finite width effects for quarkonia 
by implementing c-quark widths of 0.05 GeV 18", induc- 
ing w 0.1 GeV. The correlators show variations of the 
order of a few percent, indicating a small sensitivity to 
phenomenological (sizable) values of quarkonium widths. 
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Fig. 4. Left: cc spectral function for S— wave scattering at 
several temperatures. Right: normalized correlation function 
at several temperatures. 




Fig. 5. Same as in Fig. |l|for cc P— wave scattering. 



Fig. 6. S— wave charmonium correlators when replacing the 
"reconstructed" correlator (used for normalization) by the cal- 
culated result at 1.1 Tc. 
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